Localized end states in density modulated quantum wires and rings 
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We study finite quantum wires and rings in the presence of a charge density wave gap induced by 
a periodic modulation of the chemical potential. We show that the Tamm-Shockley bound states 
emerging at the ends of the wire are stable against weak disorder and interactions, for discrete open 
chains and for continuum systems. The low-energy physics can be mapped onto the Jackiw-Rebbi 
equations describing massive Dirac fermions and bound end states. We treat interactions via the 
continuum model and show that they increase the charge gap and further localize the end states. 
In an Aharonov-Bohm ring with weak link, the bound states give rise to an unusual 47r-peridodicity 
in the spectrum and persistent current as function of an external flux. The electrons placed in the 
two localized states on the opposite ends of the wire can interact via exchange interactions and this 
setup can be used as a double quantum dot hosting spin-qubits. 

PACS numbers: 85.35.Be, 73.63.Nm, 03.67.Lx 



Introduction. Over the last decades a number of pro- 
posals have been made for solid-state implementations 
of a quantum computer. Among these, electron spins in 
GaAs quantum dots [2 [5] are most promising candidates 
with unusually long coherence times [3]. Such dots con- 
tain typically many levels which are filled according to 
Hund's rule. Thus, the condition for a spin-qubit, which 
requires the presence of only a single unpaired electron, 
becomes challenging pi , and the scalability of such an ap- 
proach is still an open problem. In this letter we propose 
a simple setup, involving periodically modulated gates on 
top of a quantum wire (see Fig. [ij, which eventually re- 
sults in an effective double dot system. Due to the spatial 
modulation of the gate voltage the energy spectrum of 
the quantum wire acquires a charge density wave (CDW) 
gap. Recently, similarly modulated setups have been dis- 
cussed with focus on metal-insulator transitions |4] and 
transport properties in an infinite-wire superlattice [5]. 
Here, we show that the modulated quantum wire sup- 
ports localized states at each end of the wire, known as 
Tamm-Shockley bound states H] , with their energies 
lying inside the gap. These end wave-functions are well 
protected from the continuum and can host stable spin- 
qubits. 

We consider one-dimensional (ID) discrete and con- 
tinuum models and find a number of remarkable features 
for the end states resulting from the CDW modulation. 
In particular, using exact numerical diagonalization of 
the discrete open chain we analyze the stability of these 
states in the presence of a random potential and find 
that for weak disorder the end states remain stable. For 
the continuum model we consider a periodically modu- 
lated potential of the form AQCOs{kcDwx + "d), where 
Ao is the strength of the potential, kcow the CDW vec- 
tor, and 11) a constant phase. For d = 7r/2, the CDW 
phase supports zero energy bound states which are re- 
markably robust to position dependent fluctuations in 
Aq. We also show that for i? = tt/2 the model maps 




FIG. 1. The figure shows a quantum wire (black) of length 
L with negatively charged gates (blue) forming a superlattice 
potential. Due to the induced charge density modulation a 
bound state at each wire end can emerge. 

to the Jackiw-Rebbi model for massive Dirac fermions 
with midgap bound states ^6j. We treat interactions via 
fermionic and bosonic techniques and find that they pri- 
marily renormalize the gap and decrease the localization 
length. We consider end states in a ring-geometry by 
connecting them directly via tunnel junction (see Fig.|4|. 
The Aharanov-Bohm oscillation in such rings exhibits an 
unusual 47r periodicity, providing a striking signature of 
the existence of end states. Finally, we show that the 
two opposite end states serve as effective double quantum 
dot which can be used to implement quantum computing 
gates for spin-qubits. 

Lattice model. The typical lattice model for ID spin- 
less fermions in the presence of CDW modulation is de- 
scribed by ^ 

JV-l N 

H = -t"^ [c]^;^Cj+h.c.]-i-A^cos [2kcDwja+^]c^^Cj, 

3=1 3=1 

(1) 

where Cj is a fermion operator at the site j, N is the 
total number of lattice sites, t > is the hopping in- 
tegral, A > the CDW gap, kcow the CDW wave- 
vector, a the lattice constant and is an arbitrary 
phase. The energy spectrum under the constraint of 
open boundary conditions is obtained by exact numer- 
ical diagonalization and we find that the criterion for 
the existence of bound states depends on the sign of the 
potential at the beginning and end sites. For illustra- 
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FIG. 2. (a) The part of the spectra around the gap of the 
Hamiltonian given by Eq. |T|, obtained by exact diagonal- 
ization. The red bars denote two almost degenerate bound 
(midgap) states. We have chosen for the parameters t = 7, 
A = 0.8, and iV = 320. (b) One of two bound states. Plotted 
here is tp+ = ^pL + '4>R, where V'i.fl ^re states localized at the 
left (right) end of the wire. 
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FIG. 3. The dependence of the root-mean-square value 
a[Ei\ of the i-th energy level {i = 1 . . . A'^, i.e., for all en- 
ergy levels) on the standard deviation of the random disorder 
potential. 




tive purposes we have considered kcow — 7i'/4a and 
■d — 7r/2, this choice corresponds to negative potential 
at the initial two sites with the overall profile given by 
AEf=iCos[W2 + ^/2] = A(-l, 0,1, 0,-1,...). If the 
phase of the potential is chosen such that one end of the 
wire has positive whereas the other end has negative po- 
tential then only one end state is obtained. On the other 
hand, for a reflection symmetric potential profile about 
the center of a long wire (with both ends having negative 
potential), there will be two degenerate mid-gap states, 
and TpL, localized at the right and left boundaries 
resp., being the well-known Tamm-Shockley states [Hin]. 
Fig.[2|^a) shows the spectrum of an 320 site chain. Reduc- 
ing the wire length causes exponentially small splitting 
in the energies of the bound states, with the new states 
described by the symmetric and anti-symmetric combi- 
nation of "0^ and -0^. We obtain the bound states to be 
in the middle of the gap only when d = tt/2 and i ^ A. 

Disorder effects. For realistic systems, some degree 
of random disorder is unavoidable. To study this effect 
in our lattice model, we add a random on-site poten- 
tial Vic\ci. Here, Vi is taken according to a Gaussian 
distribution with zero mean and standard deviation 7. 
Fig. [3] depicts the linear dependence of the root-mean- 
square ^ (j[E.i] of the i-th energy level (i = l...A^, i.e., 
for all energy levels) on the standard deviation of the 
random disorder potential [in]. Since the slopes of the 
bound states are less than 1, we conclude that the end 
states remain gapped even for disorder strengths com- 
parable to the gap (A). As 7 is increased, Anderson 
localization sets in. We also observe as 7 is increased 
that the end states begin to mix with other (spatially) 
nearby localized states, thus effectively causing the end 
states to be more delocalized. Additionally, it is readily 
observed from Fig.[3]that the end states are more affected 
by disorder compared to all the continuum states. The 
ratio thereof depends on ^/L, since this difference is com- 
ing from the spatial localization of the end states. For a 



weak disorder, the aforementioned dependence is linear, 
while for a strong disorder the dependence becomes more 
complicated due to the emergence of Anderson localiza- 
tion. 

So far we have considered a particular realization of the 
lattice model. We next consider the continuum case, this 
limit describes the low-energy physics of a large class of 
one-dimensional lattice models with CDW (or superlat- 
tice) modulation. Recently, there has been intense activ- 
ity on exotic quantum matter, such as Majorana fermions 
(chargeless) [TTHII] and massless Weyl fermions [53] 
among others. Here, we will show that our setup allows 
for the realization of the Jackiw-Rebbi Hamiltonian 
describing a massive Dirac fermion of charge 1/2 as end 
state. 

Continuum model. We consider a quantum wire in 
the presence of a gate-induced potential with periodicity 
^CDW — I^CBW ■ For carrier densities smaller than 
the intraband energy gap only the lowest subband is oc- 
cupied. The physics of the fermion mode ^ „ {cr ^f, ]. 
is the spin index) in the lowest subband is described 
in terms of the slowly varying right TZ^{x) and left 
Ccr{x) parts and is expressed as 'i'a{x) — TZa{x)e^'"'^ + 
Ccr{x)e^^'"'^ . For an open wire, the boundary condi- 
tion 'i'cr{x = 0) = imposes the constraint [Ml [IS], 
TZa{x) = —Ca{—x). Thus, the Hamiltonian can be ex- 
pressed in terms of right movers only. 

The non-interacting Hamiltonian can be written as a 
sum of two parts, i/o — H^^'' + where the kinetic 

part can be expressed in terms of only the right moving 
fermions (the original range [0, L] now becomes [— L, L]) 
and is given by H^p = —ivp j'^j^dxTZ\,{x)dx'R-a{x) 
(summation on the spin indices is assumed) and 
the CDW term by h'^'' ~ Aq dxcos{2kcDwx + 
•d)'^\{x)'^ , with being a constant phase factor. 
Thus, i/o = (l/2)/f^da;R^-HoR<T, where R^(a;) = 
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[TZa-{x),'R-a{—x)]'^ and the Hamiltonian density Hq for 
each spin is the same and given by 

Ha = -ivpTzdx + mi{x)Tx + m2{x)Ty, (2) 

where 

mi{x) = — cos[25kx + i!?sgn(x)]Ao/2, 

m2{x) = sin[2Skx + i?sgn(a;)] Ao/2, (3) 

and Sk = kcDW — kp. li 5k = Q and the charge-density 
wave vanishes at the boundary, i.e., d = 7r/2, then it is 
easy to verify that "Hq satisfies the 'chiral symmetry' [55] 
VHf) = —HoP {V is a complex conjugation operator). 
Moreover, the eigenvalue equation (HqV'o- = ^V'ct) of the 
chiral symmetric Hq is related to the Jackiw-Rebbi equa- 
tion describing massive fermions, [6l [2^ 

-^JR^JR ^ [^^5^ _^ msgn(x)]V'''^ = er^ip'^^, (4) 

via the transformation, tp'^^ = U~^Tyi/j and H"^^ = 
U^^HaTyU, where U — exp(iT • n27r/3) and n = {i + 
j + fc)/\/3. Here, r^^y^z denote Pauli matrices act- 
ing on the spinor Tl^^x). Solving the eigenvalue equa- 
tion for L ^ vf/Aq one obtains exponentially decay- 
ing bound states tpcr ~ exp[— (Ao/2?Ji?)x] and V-'o- 
exp[— (Ao/2i;i?)(L — x)] a,t x — and x — L. Away 
from the chiral symmetry point (i9 ^ '^1'^) bound 
states still exist as long as sini? > 0, with the eigen- 
states given by V'o- ~ exp[— i(Ao exp[— j'i9]/2wi?)a:;] and 
^cr ~ exp[— j(Ao exp[— z'i?]/2wi?)(L — a;)]. For infinite 
wires the eigenvalues are degenerate and given by e = 
— Aq cos('i?)/2. However, finite length introduces overlap 
between the end states leading to an exponentially small 
splitting in the energy (see below and Fig. |4| . 

In a realistic quantum wire the gap A(x) and the phase 
•d(x) will invariably be position dependent. Assuming 
this dependence to be weak, the correction in lowest order 
in 5{x)l Aq <C 1 is given by, 

A 

Se^ d2;(5(a;)sin2i?(x)e-'^°"'"[''°]^/"^, (5) 

4wi? Jo 

where {S{x)) = and {^{x)) = {^q + S^{x)) = i?o, and 
they both vary slowly on the Fermi wavelength Xp = 
2n/kp. Thus, 6e Aq, and the bound states remain 
stable to weak perturbations. 

Interaction effects. In the following, we consider the 
effect of repulsive interactions on the end states. For sim- 
plicity, we consider spinless fermions with kp — kcDW 
and 1? = 7r/2. As usual in ID, the interactions can 
be split into forward and back scattering parts. The 
former. Hp = TTVpg4^ dx{:JiiJji: + tJiJ^:), where 
Jr = TZ'^{x)TZ{x) and Jl = C^{x)C[x), is responsible 
for the velocity renormalization [28j . v — vp{l + g^). 
On the other hand, the backscattering part, Hb = 
'^vpg2 cIx-.JlJr- at the lowest order in interaction 
renormalizes the gap. The mean-field gap A(a;) oc 



g2Vp{TZ{x)TZ^—x)) adds to the externally induced gap 
m2{x) = sgn(a;)Ao/2. We note that similar to 7712(0;), 
A(0+) = — A(0_). This can be seen by invoking the 
boundary condition, TZ{x) — —C{~x), and by express- 
ing TZ{x) = exp{i^/iTT(pfi) and C{x) = exp(— «-\/47r0i) in 
terms of the bosonic fields (I>r{x) and (phix) which them- 
selves satisfy [H], V47r(/>fl(0) = -V47r0L(O) +Tr. Thus, 
for weak interactions the bound states retain the same 
form as for the non-interacting case but with renormal- 
ized velocity and gap. To estimate the gap size we eval- 
uate the self-energy, S, using the unperturbed Green's 
function for an infinite wire, GoiiuJ,k) — (iuj — vpkr^ — 
AoTy/2)~^. In leading order, the gap renormalizes to 
(Ao/2)[l + (32/4) ln[min(A, vfL'^)/ Aq]], where A is the 
band width. Thus, the localization length, given by 
^ = {2vp/Ao){l +34 - (52/4) \n[vp/LAo]}, reduces with 
interaction. In other words, due to the repulsive inter- 
action between the continuum and the end states, the 
latter states get squeezed. 

The renormalization of the gap can be more rigor- 
ously analyzed via bosonization. Using standard proce- 
dures [29], we obtain the following form for the bosonic 
Lagrangian 

u—c,s 

+ 7r^ Vri sin[V47r(?!)n - 26kx - 

where the subscripts c, s refer to charge and spin, resp. 
The dx(f>c/s field describes the charge/spin density fluc- 
tuations and 0c/s is the conjugated field, and (p-f^i — 
{(j)c±<j)s)/V^- The Luttinger liquid parameters K^/^ and 
velocities Vc/s encode interactions, and y^^i — qAq/vf. 
The sine term denotes the coupling of up and down spin 
fermions with the external potential. As before, we as- 
sume Sk = 0. In general, there are two additional terms: 
one of them arises due to backscattering between oppo- 
site spin electrons and is given by cos{V8Tr (j)s), and the 
other, cos(4-y/7r(/)c — Akpx), describes the Umklapp scat- 
tering. However, both can be neglected as the two op- 
erators flow to zero under a renormalization group (RG) 
treatment. 

The scaling dimensions of sin(%/47r(/)-^_^), df^i = 
{Kc + Ks)/2 1, indicate that near commensurabil- 
ity [Skvp/ Aq <^ 1) the sine terms are strongly relevant. 
The parameters y^^^ have an identical flow [so as to pre- 
serve the SU{2) symmetry, this also implies Ks — 1] 
towards the strong coupling regime and yields an effec- 
tive localization length ^ ~ a(aAo/wF)^^^'^^^^'' for the 
bound state. Thus as before the role of the interactions 
is to reinforce the externally induced gap. We note that 
under RG additional terms of the type di4>cdi4>8 (where 
i — x,t) are generated, however, they are marginal and 
leave the essential physics unaltered. 

Detection. A viable approach for detecting the energy 
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The energy eigenvalues are 



FIG. 4. (a) Quantum wire (black) in an Aharonov-Bohm- 
ring geometry with negatively charged gates (blue). The 
bound states are localized on either side of the weak link 
(grey) of strength to- The energy (b) and persistent current 
j — —dF/d^ (c) dependence on the flux $/$o are plotted 
with the solid curves for the effective model (Eq. |7|) and 
with the dashed curves for the lattice model (Eq. (hj)). The 
parameters for the red and blue solid curves correspond to 
5 /to — 2.2 and (e+ + e_)/to = 6.2, and for the green curve to 
5 /to = 2.0 and (e+ + e_)/to = 6.0 (see Eq. (|8|). While for the 
dashed curves the parameters are A/to = 3.85 (red and blue) 
and A/to = 4.54 (green). The ratio A/t = 0.5 and iV = 50 is 
the same for all three dashed curves. Here, to is chosen such 
that we have a degeneracy at $/<l?o = 27r. Assuming only the 
lower bound state is filled, the persistent current shows an 
unusual 47r-periodicity as function of $/<l?o- 



splitting between the bound states is through persistent 
current measurements. For this the wire should be in a 
ring geometry so that the end states are connected to- 
gether via a tunnel junction and also large enough such 
that the energy splitting between the bound states is 
small yet the overlap of the localized wave-functions re- 
main non-zero. Such a set-up can enclose magnetic flux 
<I> inducing Aharonov-Bohm (AB) oscillations in a meso- 
scopic (phase-coherent) regime. Next consider a single 
electron placed in one of the bound states. The effective 
Hamiltonian for the spinless fermion in terms of the or- 
thogonal symmetric, |-|-), and anti-symmetric, |— ), states 
can thus be written as [7] 



Tjto COS — 



x[|-)(+|-|+>(-| 



|77)(77| + Uq sin 



$0 



(7) 



where $o = h/e is the flux quantum, e+ (e_) the en- 
ergy of the symmetric (anti-symmetric) mode, and the 
tunneling across the weak link is associated with a factor 
i^to exp(i^$/'I>o)i where v,^, — ±1 and to the tunnel- 
ing amplitude. For (anti-) clockwise tunneling we have 
/i = -|-(— ), while the sign of v depends on the relative 
sign between the wave- functions across the weak link. 
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tl/2 



± 



.y4t2 + J2_4ip^cos[$/$o]), (8) 



where 5 — \e+ — At $/<i>o = 27r the separation 

between the two eigenvalues is minimal and given by 
|2to — 5\. For large separations, the energy levels ex- 
hibit the usual 27r dependence on the flux $/$o- In 
contrast, for a flux sweep-rate w larger than |2to — 5\ 
a scenario emerges wherein an electron placed in one of 
the levels can jump to the second level and come back 
to the original one after a second 27r phase, thus exhibit- 
ing an unusual 47r-periodicity in the persistent current, 
j = —dF/d^, where F is the free energy [30]. By inde- 
pendently varying to and cu the splitting 6 can be esti- 
mated. For typical values to ~ <5 ~ 10/ieV we estimate 
j ^ O.lnA, which is of measurable size [311 I32]- For 
the observation, the phase-coherence length of the 
ring needs to exceed L. For GaAs rings, we note that 
> fiiR for sub-Kelvin temperatures [311 132] • 

The effective model, Eq. ([7|, does not take into account 
the contribution arising from the filled Fermi sea of con- 
tinuum states. However, when the number of continuum 
states below the gap is even — the states come in pairs 
with mutually canceling contributions to the current. On 
the other hand, when this number is odd, the topmost 
filled continuum state contributes to the current. Nev- 
ertheless, the amplitude of the persistent current, due to 
the end and continuum states, scale differently with the 
lattice length — the latter behaves like 1/N, while the 
former like S ^ e~^/^°. Thus, for chains with A^ ^ 1 and 
A ~ Hvp/Na, the persistent current will be dominated 
by the end states and our effective description fully ap- 
plies. The dashed curves in Fig. [3] include contributions 
from the bound states as well as the filled Fermi sea. In- 
deed we have confirmed that the contributions from the 
continuum states are two orders of magnitude less com- 
pared to those from the bound states. Finally, for the 
spinfuU case, the amplitude of j simply doubles, whereas 
the periodicity remains unchanged. 

Effective quantum dot. Similar to the discrete quan- 
tum dot states, the presence of spinful, CDW-induced, 
localized states in the quantum wire opens up an intrigu- 
ing possibility for the realization of a quantum computer 
device. These states are well separated from the con- 
tinuum and can be filled by tuning the chemical poten- 
tial to the end state level. We note that these 'quan- 
tum dots' contain automatically only one orbital level, 
and no individual gates are needed to tune them into 
a single electron regime. Due to incomplete screening 
there will be half-filling, i.e., only one state on either 
end will be filled. This is simply because once one of 
the energy levels on either end is filled, to fill the re- 
maining two levels requires additional energy to over- 
come the Coulomb repulsion. The physics of the half- 
filled state is described by the usual Hubbard model. 
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where t is the tunnehng amplitude and U is the on- 
site repulsion. For the energy hierarchy A ^ [/ t 
the effective Hamiltonian acquires the Heisenberg form, 
H = JSr ■ Sl, where J — At^/U. The effective exchange 
coupling J can be controlled by changing the gate poten- 
tial which determines the overlap between the left and 
right end modes and hence the tunneling amplitude t. 
We note that for weak overlap, t is small and U large 
making the J to be small, whereas for strong overlap 
the opposite is true [33]. By switching on and off the 
exchange constant in an appropriate sequence, the es- 
sential operations of the quantum dot, both the 'swap' 
and 'square-root-of-swap' operations can be performed, 
which, together with two single spin-qubit operations, 
enables the fundamental XOR gate [T]. 

Finite overlap between the right and the left end states 
can be ensured if their localization length ^ is on the order 
of the wire length L. This restriction yields an estimate 
for the strength of the periodically modulated external 
voltage, Aq ~ A(a/L)(^^^'=)/2, where A vp/a is the 
band width. A GaAs quantum wire with length L ~ l/im 
with approximately 10 — 20 gates requires a Fermi wave- 
length Xp ^ 50nm. And with the parameters [33], A ~ 
0.2eV, Kc = 0.8, and lattice spacing a « we obtain 
Aq ~ 0.04meV. Thus, the upper bound for temperatures 
are in the achievable range of a few hundred milli-Kelvin. 

Conclusion. We have shown that a CDW gap in a 
quantum wire can lead to bound states at the ends of the 
wire which are stable against weak disorder and interac- 
tions. They map to massive Dirac fermions desrcibed 
by the Jackiw-Rebbi model. In an AB-ring, the bound 
states lead to an unusual 47r-periodicity in the persis- 
tent current. Finally, the two opposite end states serve 
as effective double quantum dot which can be used to 
implement quantum computing gates for spin-qubits. 
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